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Abstract. We construct rational extensions of the Darboux-Poschl-Teller and isotonic po¬ 
tentials via two-step confluent Darboux transformations. The former are strictly isospectral 
to the initial potential, whereas the latter are only quasi-isospectral. Both are associated to 
new families of orthogonal polynomials, which, in the first case, depend on a continuous pa¬ 
rameter. We also prove that these extended potentials possess an enlarged shape invariance 
property. 
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1 Introduction 

Since the seminal paper of Gomez-Ullate, Kamran, and Milson [23] , which introduced the concept 
of exceptional orthogonal polynomials (EOP), the discovery of their connection with transla- 
tionally shape invariant quantum potentials (TSIP) by Quesne [41, 42], and the construction of 
infinite sets of such potentials by Odake and Sasaki [39], much progress has been made in the 
understanding of exactly solvable systems related to orthogonal polynomials (see [22] and refe¬ 
rences therein). The key tool to generate such systems is the Darboux or Darboux-Backlund 
transformation (DBT), which connects pairs of intertwined Hamiltonians. Starting from one 
primary TSIP, specific symmetries of this last select the quasi-polynomial formal eigenfunctions 
that can be used as seed functions to build chains of rationally extended potentials. The eigen¬ 
states of these extensions are then (up to a gauge factor) exceptional orthogonal polynomials, 
which, by using Crum formulas, can be expressed as Wronskians of classical orthogonal poly¬ 
nomials. The regularity properties of the chains, including degenerate chains (i.e., chains with 
repeated use of the same seed functions), are controlled by enlarged versions of the Krein-Adler 
theorem [3, 17, 22, 32, 47]. For some chains, the extended potentials share the same shape inva¬ 
riance properties as the primary potential [24, 25, 42, 43]. With other choices of seed functions, 
the resulting potentials possess an enlarged shape invariance property [26, 44, 45]. 

Until now, the chains of extensions were “rigid” in the sense that they were uniquely de¬ 
termined by the tuple of associated seed functions. Very recently, with B. Bagchi [4], we ob¬ 
tained new rational extensions of the Darboux-Poschl-Teller potential (based on the so-called 
para-Jacobi polynomials [9]), which depend on a free parameter and can then be modulated 

*This paper is a contribution to the Special Issue on Exact Solvability and Symmetry Avatars in honour of 
Luc Vinet. The full collection is available at http://www.emis.de/journals/SIGMA/ESSA2014.html 
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continuously. The eigenstates of these extended potentials are associated to new families of or¬ 
thogonal polynomials that are, in a broad sense, exceptional para-Jacobi polynomials and which 
depend on a free continuous parameter. 

In this paper, we consider the possibility of building new rational extensions of two confining 
TSIP, namely the trigonometric Darboux-Poschl-Teller (TDPT) and isotonic potentials, via 
confluent chains of DBT, that is chains of DBT in which the spectral parameters of the different 
seed functions converge to the same value. It has to be noticed that it is precisely by using 
such confluent chains applied to the constant potential and considering the associated rational 
extensions that Adler and Moser built the Burchnall-Chaundy polynomials [8] in their seminal 
paper on the rational solutions of the KdV equation [2] . 

The possibility of considering two successive factorization energies tending towards a common 
real value in a chain of Darboux transformations was first considered in the framework of phase- 
equivalent potential construction [5], then extended to a class of potentials defined on the 
line [49]. These approaches generalized to arbitrary bound-state energies a procedure already 
known for the ground state and the first few excited states [13, 31]. An independent proposal 
was made wherein the terminology “confluent” algorithm was introduced and some partners 
of the free particle and the harmonic oscillator were exhibited [37]. The confluent algorithm 
was then studied in more general terms and applied to the free particle, one-soliton well, and 
harmonic oscillator [19]. It was also considered in a general construction of all possible first- and 
second-order partners of the TDPT potential [12]. The “hyperconfluent” third-order algorithm, 
wherein the three factorization energies converge to the same value, was analyzed [21]. Some 
Wronskian formulas applicable to the confluent case were also derived [7, 20, 48]. 

The present work differs from the previous ones devoted to confluent chains by the restriction 
to final potentials that are rational extensions of the initial ones. As a consequence of this 
condition, some parameters of the latter may have to be chosen integer. 

After recalling in Section 2 the basic elements concerning the Darboux-Backlund transfor¬ 
mations, we review in Section 3 the concept of confluent chains of DBT. We show in particular 
that the confluent chains of arbitrary order can be generated within the standard frame of 
(completed) DBT chains, giving rise to multiparameter dependent extensions. 

In Section 4, applying two-step confluent chains of DBT for which the seed functions are 
eigenstates, we build regular rational extensions of the TDPT with appropriate parameters. 
The extended potentials depend on a continuous parameter and are strictly isospectral to the 
initial potential. The eigenstates form new families of orthogonal polynomials, which have a free 
parameter dependence. We exhibit particular examples and prove that the extended potentials 
present an enlarged shape invariance property, in which the parameter transformation acts in 
a nontrivial way on the supplementary parameter. 

In Section 5, we make the same construction for the isotonic system. In contrast with the 
TDPT case, the regular rational extensions do not depend on any supplementary degree of 
freedom and we only have quasi-isospectrality between the extended potentials and the original 
one. We also furnish explicit examples of extensions and establish their enlarged shape invariance 
property. Section 6 contains some final comments. 

2 Darboux—Backlund transformations: basic elements 

We consider a one-dimensional Hamiltonian H = —d 2 /dx 2 - t-H(x), x G I C M, and the associated 
Schrodinger equation 

^(*) + (E\ - V(x))^{x) = 0, (2.1) 

ip\(x) being a formal eigenfunction of H for the eigenvalue E\. In the following, we suppose 
that, with Dirichlet boundary conditions on /, H admits a discrete spectrum of energies and 
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eigenstates (E n , V ; n)ne{o,...,n max }CN) where, without loss of generality, we can always suppose that 
the ground level of H is at zero (Eq = 0). 

The Riccati-Schrodinger (RS) function w\(x) = —ip' x (x)/ip\(x) associated to ipx satisfies the 
corresponding Riccati-Schrodinger equation [28] 

-w'\{x) + w\(x) = V(x) - E x - ( 2 . 2 ) 

From any solution ip v (or equivalently w v ), we can build a Darboux-Backlund transformation 
(DBT) A(w u ) defined as [10, 11, 15, 16, 28] 

w\(x) A{ -^ v) w^\x) = -w v (x) + (E x - E v )/(w v (x) - w x (x)), 

ip\{x) A H l ' ) ipx\x) = exp J dxw^\x) \ ~ A(w„)ip\(x), A ^ v, (2.3) 

where A(w u ) is a first-order differential operator given by 
A(w v ) = d/dx + w v (x). 

ip^ and are respectively solutions of the Schrodinger and RS equations with the same 
energy Ex as in equations (2.1) and (2.2), but with a modified potential 

V^\x) = V(x) + 2 w' u (x), (2.4) 

that we call an extension of V(x). For the associated Hamiltonian = —d 2 /dx 2 + V^ v \x), 
we have the factorizations 

H {u) = A(w u )A+(w u ) + E V , H = 4+K)%) + E v , 


with 


ipx(x) ~ A + (w„)ip^\x). 

The function ip^ in equation (2.3) can then be rewritten as the Darboux-Crum formula 

(2.5) 

where W(y\,... ,y m \x) denotes the Wronskian of the family of functions y±,... ,y m , 


,{v), , W(ip v ,ip x \x) 

« (I) ~ M*) 


W(yi, ...,y m \x) 


yi(x) ... y m ( x) 

(m—l) / x (m—1) / x 

y\ (x) ... ym \x) 


The eigenfunction ip v is called the seed function of the DBT A(w u ) and and ip^ are 
the Darboux transforms of V and ipx, respectively. 

Note that A(w u ) annihilates ip v and, consequently, equations (2.3) and (2.5) allow to obtain 
an eigenfunction of V ^ for the eigenvalue Ex only when A 7 - u. Nevertheless, we can readily 
verify that 1 /ip u (x) is such an eigenfunction. By extension, we then define the “image” by A(w u ) 
of the seed eigenfunction ip,, itself as 

^v\x) ~ 


ip v (x)' 


( 2 . 6 ) 
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At the formal level, the DBT can be straightforwardly iterated and a chain of m DBT can 
be simply described by the following scheme 

A (^l) ,( Vl ) A { W "2 l) ) (JV 2 ) A i. W "rrT J) ) , N) 

V’A 5 * ~ Wx ■" w > 


y ^ yOD v ^4 ' yCYd 




y(JVm) ^ 


where IVj denotes the j-uple (iq,..., Vj) (with A r i = zq), which completely characterizes the 
chain. We denote by (N m , v m +ii • ■ •, zqn+fc) the chain obtained by adding to the chain N m the 
DBT associated to the successive eigenfunctions n ■ • • > 

- g an e ig en f unc ti on associated to the eigenvalue Ax of the potential (see equation (2.4)) 

m 

V {Nm \x) = V(x) + 2^2(wl^ j - 1 \x))' (2.7) 

7=1 


and can be written as (cf. equations (2.3) and (2.5)) 

= A{ W ^-^)^- l \x) = A^w^-A) • ■■A(w Vl )^x{x). 


( 2 . 8 ) 


A chain is non-degenerate if all the spectral indices zq of the chain N m are distinct and is 
degenerate if some of them are repeated in the chain. For non-degenerate chains, Crum has 
derived very useful formulas for the extended potentials and their eigenfunctions in terms of 
Wronskians of eigenfunctions of the initial potential [14]. 

Crum’s formulas. When all the vj and A are distinct, we have 


^N m ) 


(x) 


W( JVm ’ A )(x) 
W( Nm \x) 


(2.9) 


and 


V {N ™\x) = V(x) - 2(log W iN ™\x))'', 
where W^ Nrn "> (x) = W (z/Vi,..., i/} Vm \ x). 


( 2 . 10 ) 


3 Confluent chains of DBT 

The single-confluent limit of a chain of DBT N m is obtained when all the spectral indices in¬ 
tend simultaneously to the same value Vj -A- u, V j e {1,... ,m} (in the following, we consider 
only single-confluent chains and than omit the adjective “single”). 

3.1 Two-step confluent chains 

We consider a chain of two DBT = (zq, v 2 ), which, in the non-degenerate case zq 7 ^ zq, gives 
(see equations (2.7), (2.8), (2.9), and (2.10)) 

y(^)( x ) = V(x) - 2(logfC , ' 1 ’ I ' 2 )(x)) // = V(x) -2 [{E V2 - E vl )/{w V2 {x) - w Ul (x))]', 
^Xx) = {Wui(x) - w U2 (x))ip U2 (x) = W(il) ul ,4> V2 I x)/ip ui (x). (3.1) 

Note that in the degenerate case zq = zq, we have (see equation (2.6)) 

^i 1 X x ) = l/^ Vl (x), w^\x) = -w Ul (x), 




Confluent Chains of DBT: Enlarged Shape Invariance and New Orthogonal Polynomials 


5 


and by applying the DBT A(w = A(—w Ul ) 
V^ vi \x) = V^\x) - 2w' vi (x) = V(x). 


to V (vi) , we recover simply the initial potential V, 


The confluent case corresponds to the limit ^2 —> v\. As proven by Fernandez et al. [7, 20], 
the confluent extended potential and its eigenstates admit the following integral representations 


V^\x) = V(x ) - 2 (log ( / dtifi.it) - Wo 


’ XQ 


and 


^ 1,Vl \x) = ( E V1 - E k )ip k (x) - 


Vd (®) 




Ix, j^(t)-Wo 


(3.2) 


(3.3) 


Both depend on an arbitrary real parameter Wo and for an adapted range of Wo values, the 
extended potential is regular. In fact, the formula for the potential (3.2) already appears in 
many previous works, for instance in a 1986 paper of Luban and Pursey [33] and a few years 
later in [31]. The Matveev formulas [27, 34, 35] for the two-step case 


V^ Vl \ x ) = V(x ) - 2 


log W 



/ dil>v{x) 

V dE v 


V=l/\ 



= W 



/ d^ v {x) 

V dE„ 




( d'ipv(x) 

V dE u 



(3.4) 


which express the confluent extension and its eigenstates in terms of generalized Wronskians 
(which are in fact two-way, or double Wronskians [51]) can be viewed as associated to a particular 
choice of the Wo constant. Indeed, if we consider the indexed family of RS functions w v (x) as 
satisfying a prescribed initial condition in xo, we have, in the confluent limit, w U2 (x) —>■ w vi (x). 
It results from equation (3.1) that 


W—> W = r W -2(i/(^) 

But we can readily verify that 


/ du>v(x) 

V dE u 




log W 



f d'tj) v {x) 

V dE„ 



/ 


and since in this case we also have [36] 


dw u (x)\ 

dK ) 


1 

(®) 





we see that equation (3.4) corresponds to equations (3.2) and (3.3) with Wo = 0. 

It has to be noticed that the degenerate extension V'- Vl,I/1 >{x) can be recovered from the 
confluent one V^ ui,ui \x) by taking the singular limit value Wo —> 00 . 

Fernandez et al. used these formulas to generate new second-order SUSY partners of the free 
particle, the Kepler-Coulomb, and the single-gap Lame potentials [7, 20]. 

The preceding results can be in fact integrated within the standard DBT scheme simply 
using the DBT in its completed form (see equation (2.6)) as in [13, 31]. Indeed, by applying the 
DBT A(w v ), we generate first the one-step (possibly singular) extension 

V^ v \x) = V(x) + 2 w' v (x) = V(x) - 2(log ij) v {x))". 
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Since ip^ = 1/Vv is an eigenfunction of for the eigenvalue E u , the most general eigen¬ 
function (up to a multiplicative factor) of for the same eigenvalue is 


^\x-X 1 ) = ^\x) A 1 + 


dt- 


A i + fx o 


Ixo (V’^(i )) 2 
and the corresponding RS function is 

% + 1! dtyKt) 




Ai E 


(3.5) 


wM(s; Ar) = - 


log 


1pv(x) 


= —w v (x) - 


tpl(x) 


A i + Sxo dtyl(t)' 


We can now use this general solution as seed function for the second DBT. Then ap¬ 
plying A(W^ y ) to V^ u \ we obtain the following second extension 


V^\x- Ar) = V^\x) - 2(log^(x; Ar))" = V(x) - 2[log (^(s)^M(x; A,))]' 


= F(x) - 2 


log (Ai + f dt'ipl(t) 

\ Jx 0 


= F(x) - 2 


ipl(x) 


X l+fx 0 dt 'ft( t ) t 


(3.6) 


and we recover the first Fernandez formula (3.2) with Ai = —Wq. As for the eigenfunctions 
of , they are given by (fj, / v) 


^\x-\i) = A(W^)^\x) = (E u - E^(x) 


W (ipy, 1 X) 
^(sjAi) 


(3.7) 


3.2 General multi-step confluent chains 

Fernandez and Salinas-Hernandez [21] have also considered the so-called “hyperconfluent” case 
corresponding to a three-step confluent DBT, for which they have extended the previous for¬ 
mulas, equations (3.2) and (3.3). For these three-step extensions, the potential depends on two 
arbitrary real parameters. 

In fact, the preceding analysis allows to obtain integral formulas “a la Fernandez” for chains 
of arbitrary order in a very simple way. In the following, the symbol (is 1 ) means {y,... ,v). 

'-V-' 

l times 

In the three-step case, the image of (see equation (3.5)) by the DBT A(W^) being 
^ = l/’I'i^, the general eigenfunction of V^ 2 \x\ Ai) (see equation (3.6)) associated to the 
eigenvalue E u is 


*fXx;A 2 ) = 


A 2 + /; o dt(^ ) (t;A 1 )) 
^(xjAi) 


Ai, A 2 E M, 


where we have used the notation A m = (Ai,..., A m ). The next extension generated by the DBT 
A(wl u ■*) is then 


F^(x;A 2 ) = R(x)-2[log (^(s)®M(s; \i)*¥\x- A 2 ))]" 

= V(x) - 2(\og'ip v (x))" - 2 (\ 2 + f dt(^l u \t] Ai)) 2 

V dx 0 

We recover the “hyperconfluent” third-order superpartner of V(x) as obtained by Fernandez 
and Salinas-Hernandez [21]. Within this scheme, the generalization is immediate and repea¬ 
ting the procedure m times, we obtain for the hyper confluent m th -order extension of V{x) the 
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expression 


m— 1 


771—1 


A m _r) = F(x) - 2 log J] T^)(x; A,) = P(x) - 2 £ (log ^(x; A,))", 

\ 3=0 J 3=0 

with the following recurrence relation for the successive seed functions 

a k + r A fc _o ) 2 


^\x;A k ) = 


x 0 






(x; A fc _x) 


where Tj^x; Ao) = Vv(x). 

In other words, whenever m is even (m = 2k) 


F (i/ 2 L) (x; A 2 fe _i) = V(x) - 2 |log 
and whenever m is odd (m = 2k + 1 ) 
V {u2k+1 \x-A 2 k) = V(x) - 2 {log 


7c-1 


U ( A 2 z+i + / \t]A 2 i))' 

U=o v 1/3:0 


iWn ^f / J) (*;a 2 z_i))' 

z=i ^ 1/3:0 


The eigenstates of V^ m ' ) can be obtained by successive applications of the A{wl v ^) operators 
(the product being ordered in decreasing order), 

m 

\ x 'i A m-l) = n ^■{ W u l/l 1) )'^k(x). 

1=1 

A direct application of the Crum Wronskian formulas [14] to these general, parameter- 
dependent, confluent extensions is obviously not possible and, as mentioned above, the Matveev 
formulas [27, 34, 35] correspond only to a particular choice of the A j parameters. Nevertheless, 
these extended potentials are amenable to other (standard) Wronskian formulas [7, 20, 48]. 

In the following, we limit our analysis to the case of two-step DBT. We are interested in 
the possibility of building regular and rational extensions with such confluent chains, which 
turns out to be possible for the trigonometric Darboux-Poschl-Teller (TDPT) potential and the 
isotonic potential. 


4 Two-step confluent rational extensions of the trigonometric 
Darboux—Poschl— 1 Teller (TDPT) potential 

4.1 General scheme 

The trigonometric Darboux-Poschl-Teller (TDPT) potential (with zero ground-state energy) is 
defined on x e ]0, vr /2 [ by 

T „ ^ _ (« + l/2)(a-l/2) , (/3 + l/2)(/3 - 1/2) , , , in2 

V (x, a, p) — -;— y -1- y - [ot + p + 1) , 

sin x — ~ 


COS 2 X 


with a, j3 > 1 / 2 . 
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Its physical spectrum, associated to the asymptotic Dirichlet boundary conditions 

^(0 + ;a,/3) = 0 = V' 

is given in terms of Jacobi polynomials [18, 29, 50] 

p(a,P)( \ (-l) w r(n + f3 + 1) v - ' / sk(n\ r(n + a + f3 + 1 + k) k 

n n!r(n + a + /3 + l) J VV 2 fc r(/? +1 + k) 1 

by 

E n (a, (3) = ( a n + f3 n + l ) 2 — (a + j3 + l ) 2 = 4n(a + j3 + 1 + n), 
i/> n (x;a,/3) = 'it>o(x;a,P)p( a ’P\z), n E N, 

with z = cos2x E ]—1,1[, (c^n, /3 n ) = (a + n, (3 + n ), and 

^o(x; a, /?) = (! — ^)(«+i/2)/2 (1 + z )(P+mi\ 


In the following, in order to get rational extensions, we consider the case where a and (3 are 
integers: a = N>l, (3 = M>1. 

If we choose as initial seed function an eigenstate ip n (x\N,M) of V(x;N,M), by taking 
xq = 7t/2 (zq = —1), the quantity 


Qi NM \z)= r dtrt&N,M) = ~ i dc(i-o N (i+o M (p^ M \oy 

J 7T/2 


1 

2 .i-i 


is a polynomial of degree IV + M + 2n + 1 in z with [18, 50] 

= f 1 


Qr 


> (1) = - \ 1 1 dC(i - 0*(1 + C) M {P^ N ’ M) (0) : 

2 iV+M (n + N)\(n + M)\ 


2n + N + M + l n\(n + N + M)\ 
Note the following recurrence 

4n 


QS 1 ,M+ 1 ) (i) = 




n + N + M + 1 

From equation (3.6), we then obtain for the confluent two-step extension V^ n2 \ 
V^\x;N,M, Ar) = V( X] N,M) - 2[log (A, + Qf’ M) (*))]" 

1 ' ' <n A 1 + (?r- M| ( 4 ) 


(4.1) 


(4.2) 


where V^ n2 \x] IV, M, Ai) constitutes a rational extension of V(x) (in the z variable). Q^’ M \z) 
is strictly decreasing on the interval ] — 1,1 [ and keeps a negative value, varying from 0 to 
Qn*’M)( 1) < 0 when z runs through ] —1,1[. Consequently, when 

A 1 E]-oo,0]u]-Qi 7 V ’ M )(l),+oo[, (4.3) 


then Ai + Q^’ M \z) keeps a constant sign, strictly negative or strictly positive respectively, and 
V^ n2 \x] N , M, Ai) is regular. 
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In this case, its eigenfunctions for k ^ n are given by (see equations (3.5) and (3.7)) 


N, M, Ar) = [E n (N, M ) - E k (N, M)]^( x) - 


with 


Wjlpn^k 1 x) 
*&°(®;IV,M, Ap 


*<">(*; N, M, A,) = (1 - ,)-(w/2)/2 (1 + zr (M+l/2)/2 ^ + g' M> (*) 


PP M, ( 2 ) 


and 


IE(^n, ^ I x) = -(1 - ^^(1 + ^) M+1 PiJ M) (^), 

where (Pi“’ /3 ^(z) = 0) 


P 


N k M \z ) = (fc + IV + M + l)Pf M )(z)p£ +1 ' M+1) ( Z ) 

— (n + IV + M + i)p^+ 1 ^+ 1 ) ( z ) P WM ) ( 2 ) 


is an exceptional Jacobi polynomial in the broad sense of the term. 

Hence, for k ^ to, 

IV, M, Ai) = [E n (N, M ) - £ fc (IV, M)]^ fc (s; N, M ) 

#,M), x p (JV,M), x 

+ (1 - ^)(3JV+5/2)/2^ + z ^3M+5/2)/2 n,k n 

Al + Qn ’ (2) 


that is 


Pif.ljz-, A,) 


if\x-,N,M,\P) = (1 - 2)< N+1 «/2 (1 + 2 )(*W/2)/2 , 

Ai + Qn ’ ; (z) 


where 


AJStte Ai) = 4(n -k)(n + k + N + M + 1)P< W ' M, ( 2 )[A, + 

+(i - 2 ) w+i (i+rtfwpr’w. 


Moreover 

„2 


4 n2 )(s;IV,M, Ar) = lM n )(ar;IV,M, Ar) 

= (1 _ 2 )( 7 V+ 1 / 2 )/ 2 (l + z )(M+ 1 / 2)/2 


p(JV,M) 

n 


(z) 


^+Q ( n NM \z) 

which is a normalizable eigenstate of I/( n2 )(x; IV, M, Ai). This corresponds to defining 

Pnm^^) = P^ M \z)- 

; (x\ IV, M, Ai) tends to zero at z = — 1 and z = 1 (i.e., x = tt /2 and x = 0) and is then 
an admissible eigenstate of V^ n ~\x\ N, M, Ai) for every k > 0. The potentials V(x; N, M) and 
y(n )(aj; IV, M, Ai) are therefore strictly isospectral. 
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/ 2 \ 

The orthogonality conditions between eigenstates imply that the Ai) constitute 

a family of orthogonal polynomials (indexed by k 6 N) on ]—1,1[ with respect to the measure 


( n 2 ) / \ \ 

M N,M\ Z ' A 1 ) 


1 (1 - z) N {l + z) M 


It is worth observing here that the confluent two-step extension V^ n2 ^(x; N, M, Ai) may be 
considered as a special case of one of those with general parameters that have been built by 
Contreras-Astorga and Fernandez in [12, Section 3.2.3(b)] (namely the third one given in equa¬ 
tion (3.43)), whenever their parameters A, v assume the half-integer values A = N + 1/2 and 
v = M + 1/2. 


4.2 Examples n = 0 

We consider the n = 0 case. Then 


Qi N ’ M) 


(*) 


-\ f da i-o n (i+o m 


N 

~(z+ i) m+i J2 

k =0 


2 J V -fe-l(_ 1 )fe /jyX 

M + k + 1 \k) 


(z + l) k 


and 


V ( ° 2 \x-N,M, Ai) 


0 (1 - z) 2N+1 ( 1 + z) 2M+1 


Moreover 


if’(*1 N, M, A,) = (1 - 2 )< K+1 /2)/2 (1 + ,)(« + l/2)/2 , 




where 


P 


NM, k ^ a i) =+ n + m + i)^f’ m) w(Ai +Qr ,M) (^)) 


+ (1 - ^^(l + z) M+1 P^’ M \z), k / 0, 

p(° 2 ) \ \ _ 1 

with (P^[’ M \z) = 0) 

P$ M \z) = (k + N + M + 1 )P^ 1,M+1 \z). 

4.2.1 The N = M = 1 case 

The results read 


q£’%) — --(^ + 1) 2 _ ^ 

y(° 2 \x; 1,1, Ai) = V(x; 1,1) — 12 


<1 - - 2 ) 


Ai - |(z + l) 2 + Uz + l) 3 
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+ 2 


(1 - z 2 ) 


2\3 


[Ai - \{z+ l ) 2 + \(z + l) 3 ]' 


or 


P(° 2 )(x; 1,1, Ax) = + 3 2 - 9 - 12 

4 sin x 4 cos z x 

sin 6 2x 

+ 2 


sin 2 2x cos 2x 


Ai — 2 cos 4 x + | cos 6 x 


(Ai — 2 cos 4 x + | cos 6 x )' 


with Ai E ]—oo, 0] U ]1, +oo[, and 


) (s;l,l,Ai) = (l-z) 3/4 (l + *) 3/4 -- T 


d '2 (-*;>■) 


Ai — | {z + l ) 2 + l(z + l ) 3 


where 


Ai - + l ) 2 + g( z + !) 3 


P i 0 ii(^ A i) = -4^ + 3) J P fe (1 ’ 1) W 

+ (A: + 3)(1 — z) 2 (l + z) 2 P^^(z), k / 0, 
P 1 ( S(z;A 1 ) = l. 


4.2.2 The AT = 2, M = 1 case 

The results read 

= — (^ + l) 2 _ + 1) + g(^ + l) 2 ^ , 

C<» 2 +; 2, 1, A,) = V(z; 2,1) - 4 (1 + 4 *' )(1 “ z) ' (l + z) 

+ 2 


or 


Ai-(z + l ) 2 + §(z + l) 3 -l(z + l ) 4 
(l-z) 5 (l + z ) 3 


[Ai - (z + l ) 2 + §(z + l ) 3 - |(z + l) 4 ] 


E(° 2 )(x; 2,1, Ai) = + . 3 . - 16 - 32 


4 sin 2 x 4 cos 2 x 

+ 512 


(1 — 4 cos 2x) sin 4 x cos 2 x 
Ai — 4 cos 4 x + ¥ cos 6 x — 2 cos 8 


x 


sin 10 x cos 6 x 


(Ai — 4 cos 4 x + ^ cos 6 x — 2 cos 8 x )" 


with Ai 6 ]—oo, 0] U ]|, +oo[, and 


+ ,) (.t;2,1,A.) = (1-+ /4 (1++ /4 




Ai — (z + l) 2 + §(z + l) 3 _ §(+ + l) 4 


where 


-^ 2,1 ,l( z ’ Ai) — —4fc(/c + 4)P^ ’ •’(z) ^Ai — (z + 1 )" + -(z + l ) 3 — -(z + l) 4 ^ 

+ (k + 4)(1 - z) 3 (l + z) 2 P fe (3 f (z), k + 0, 

P 2 (0 1 >z;A 1 ) = 1. 
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4.2.3 The N = 1, M = 2 case 

The results read 

Q?’ 2) (*) = _ ( z + 1 ) 3 ^3 “ g( z + 1)^ > 


+ 2 


(1 -z) 3 {l + z) 5 


[Ai — \{z + l) 3 + l(z + l) 4 ] 


or 


V^° \x; 1, 2, Ai) = - r ~2 h 


3 15 _ „„ (1 — 4 cos 2x) sin 2 x cos 4 x 


4 sin 2 x 4 cos 2 x 
+ 512 


- 16 + 32- 


Ai — | cos 6 x + 2 cos 8 ; 


sin 6 x cos 10 x 


(Ai — | cos 6 x + 2 cos 8 x)' 


with Ai £ ]— oo, 0] U ]|, +oo[, and 


^f>(x;l,2,Ai) = (l-+ /1 (l++ /4 




Ai - Uz + l ) 3 + jUz + l) 4 ’ 


where 


-Pi,2,l( z > ^i) _ ~4k(k + 4 : )p£ ' \z) ^Ai — ~(z + l) 3 + -(z + 1 ) ^ 
+ (k + 4)(1 - z) 2 (l + z) 3 pj£?(z), k + 0, 
^i (0 2,o(^A 1 ) = 1. 


4.3 Shape invariance of the two-step confluent rational extensions 
of the TDPT potential 

U(x;IV, Af) is a translationally shape invariant potential with a SUSY partner 
U S usy(+ N , M ) = V W (x; N , Af) = U(x; N + 1, M + 1) + E ^ N , Af) 
and (the coefficient is readily established by a direct calculation) 


4°)(x;1V,M) 


Wj^Q^n [ x) _ / 


En ( N , M ) \ 

4 n J 


Ipn— l(+ N + 1, Af + 1). 


Since V^ n2 \x; N, Af, Ai) and V(x; N, Af) are strictly isospectral, it is natural to wonder 
whether l+ n “)(x; IV, Af, Ai) shares the same invariance property as V (x; N, Af). 

Consider the SUSY partner of l+ n )(x; IV, Af, Ai). It is given by 

^susy(*; n .Ai) = U (n2) (x; IV, Af, A,) - 2( log vf 2) (x; IV, Af, A,))" 

= U(x; IV, Af) - 2 [ log ty„(x; N , Af)W( X ; IV, Af, Ai)$ n2) (x; IV, Af, Ai))]", 

where, for n > 1, 


4 n2) (*; ^ A i) = £«(W, N , Af) 


Wilpn^O I x) 
4i n) (x;JV,Af,Ai) 
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and, for n = 0, 

$° 2 ) (x;N,M, Ai) = 


^ 0 ) (x;N,M, Ar)' 


Consequently, in the n = 0 case, we have 

^susyOc N, M, Ai) = V(x- N, M) 

- 2 [log (Mx\N, M) < 0 ) (s; IV, M, Ai)^(®; IV, M, Ai))]' 
= P(x; IV, M) - 2[log V , o(^5 IV, M)]", 

that is 

^susyOc iV, M, Ai) = IV, M) = V(x; N + 1,M + 1) + Ei(IV, M). 

Furthermore, in the n > 1 case, we get 

^susy(*; ^ M . A i) = IV, M ) - 2[log ^(®; N, M)}" 

- 2 | log 

= V(x; N + 1, M + 1) + Ei (IV, M) 

~ 2 (log 


E n (IV, (*; IV, M, Ai) + ^°’^ ra ^ ) ) ) V’nfo IV, M) 


E n (N, M)^\x- N , M, Ai) - IV + 1, M + 1) 

4 n 

+ log[^n(®;E',M)]| . 

More precisely 

^SUSY^; N > M , A i) = ^(®; IV + 1, M + 1) - 2[log^ n _i(x; IV + 1, M + 1)]" + Ed (IV, M) 


- 2 log 


Aj + Q { 'n' M \z ) — il> n {x] IV, IV + 1, M + l)/4n \ 

V’ n -i(^;iv +i,M + 1 ) y 

= F (n_ 1 ) (x;IV + l,Af + 1) + E 1 (N,M) 

_ 2 A Aj + Qn Ar ’ M ^) — IV, M)ifj n —i{x'i IV + 1, M + l)/4n \ 

° S ^ n _i(s;IV+l,M+l) J 

If, for an appropriate constant C, the following condition 

Ai + - ^:W)A-ifeJV + i.M + i) = C(A , + qS ..m +1 , (2)) (4 4) 

is satisfied, then 

Ai + QL N,M> (z)- v/>„(x;N,M)i/> n - 1 (x;N+ 1,M + l)/4n ,/, 

-^_ l( , ;JV+1 . M+1) -= C *-' feJV + l.M + l.A,) 

and we obtain an enlarged shape invariance property 

V^ Y {x-,N,M,\ 1 ) = V« n -V 2 \x-,N+l,M + l,\' 1 ) + E 1 (N,M). (4.5) 

The preceding condition (4.4) can be rewritten as 

A(z)-B(z) = CX' 1 -X 1 , Vz €]-!,![ (4.6) 
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with 


A(z) = Q^’ m \z) - CQ^_X lM+1 \z 
_ 1 
“ 2 


d([c(i - 0^(1 + c) m+ 1 (pS 

-(i-() s (i + 0 M (pf’ tf) (0) 


(iV+l,M+l)^2 


and 


B(z) = ipnix] N,M)i/j n _i(x-,N + 1, M + l)/4n 

= (1 ~ z) N+1 ( 1 + ^) M+1 p(Af,M)^^p(^+l,M+l)^ 

Equation (4.6) is equivalent to 

d . d . . 

sY (z) = * B(z) - 


where 4zA(z) is given by 


-A(z) = C(l-z 2 )(Pi N t 1 ’ M+1 \z)) 2 -(P} l N ’ M \z)y 


(1 - z) N (1 + z) M dz 

As for 2 feB(z), using the derivation formula [18, 50] 

d p (N,M)(~\ - N + M + n + 1 p (JV+l,M+l), v 
dz n [) ~ 2 - 1 [Zh 

it can be expressed as 

_ ^ _ —B(z) ~ (l- z 2 )p( N ’ M Uz)—P {N+1 ’ M+1) (z) 

(1 -zr(l + z)M d z [Z) ~ [1 Z)Fn dz n ~ x {) 

+ (1 - z 2 ) W + M 2 + ” + 1 (p!” +1,M+1> (z)f 

+ [(M — N) — z(N + M + 2)]Pf’ M) (z)pf 1 +1 ’ M+1) ( Z ), 


that is 


4n 


(1 - 2 )^( 1 + 


«W = (i - ^ ) " + M + " + i (P (ywi) W )^ 


+ [(M — N) — z(N + M + 2 )}P( N ’ M \z)P^ 1 ’ M+1 \z). 

But the differential equation satisfied by the Jacobi polynomials is [18, 50] 

(1 - z 2 ) ^P^ N ’ M) {z) + m -N)~ z(N + M + 2)}j-P^ N ’ M \z) 

= —n(N + M + n + 1 )P { n N ' M \z), 
which, combined with equation (4.8), gives 

(1 - z 2 ) ^P^! + 1 ’ M+ 1 ) (z) + [(M — N) — z(N + M + 2)\pj£+ 1>M+1 \z) = -2 nPf M >( 2 ). 


(4.7) 


(4.8) 
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Consequently 

4n 


(1 -z) N (l + z) M dz 


B(z) 


= -2 n(P+"\z)f + (1 - , 2 )« + M_fn+ i(p ,« + 1 ,M + 1 , ( 2 )) 2 i 


or 


d 


(1 -z) N (l + z) M dz 


B(z) 


= (i - 0 JV+A 4,t n+ 1 P"-f wl)(z) > 2 - 

To satisfy the equality (4.7), we then must choose 
N + M + n + 1 


C = 


4 n 


In this case, equation (4.6) simply becomes 

A(z) - B(z) = -—-- Ai. 

In the limit z —> 1~, we obviously get 

B(l) = 0 , 

and (see equations (4.1) and (4.2)) 

21(1) = QW">(1) - JV + M J" + 1 Q 1-T M+1) (1) = 0. 

which implies 


4n 


Al N + M + n + 1^ 1 ' 


(4.9) 


We conclude that for n > 1, the two-step confluent rational extensions of the TDPT potential 
satisfy the enlarged shape invariance property (4.5) with A) given in equation (4.9). Note that 
the latter relation ensures that the domain of Ai values for which V^ n2 \x; IV, M, Ai) is regular 
corresponds exactly to the domain of A) values for which I/(( n_ 1 h)(x; N + 1, M + 1, A)) is also 
regular (see equations (4.2) and (4.3)). 

In the n = N = M = 1 case, such a property can be directly verified as follows. We get 


P 


( 2 , 2 ) 


(*) = 1, P[ 1 ’ 1 \z) = 2z, 


C = \ , as well as 

A(z) = Q^\z) - Q^’ 2 \z) = ^ dC((l - C 2 ) 2 - 4C 2 (1 - C 2 )) 
= l - £ d({ 1 - 6C 2 + 5C 4 ) = \z{ 1 - z 2 ) 2 


and 


(l-z 2 ) 2 

fl(z) = A-—2a;. 

4 

Hence the identity 

A(z) = B(z), Vz€]~ 1,1[, 

is satisfied, which implies that 

^susy(*; !» A i) = 2,2, Ai) + Ek(\, 1). 
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5 Two-step confluent rational extensions 
of the isotonic potential 

5.1 General scheme 

The isotonic oscillator potential (with zero ground-state energy E$ = 0) is defined on the positive 
half-line ]0, +oo[ by 


V (x; uj, a) = —x^ + 


uj 2 2 (a + 1/2) (a — 1/2) 


— Lo(a + 1), |a| > 1/2. 


If we add Dirichlet boundary conditions at zero and infinity and if we suppose a > 1/2, it 
has the following spectrum (z = ujx 2 /2) 


E n (uj) = 2nuo, ^ n (x; uj, a) = z^ a+1 ^ 2 ^ 2 e z / 2 L™{z), n> 0, 


where 


KO) = Y 


(—l) fc (n + a) ■ ■ ■ (k + 1 + a) 


k =o 


k\ 


(n — k)\ 


(5.1) 


is the usual Laguerre polynomial [18, 50]. 

In the following, in order to get rational extensions, we consider the case where a is an integer: 
a = N > 1. 

If we choose as initial seed function an (unnormalized) eigenstate uj , N) of V(x;uj, N), 
by taking xq = 0 we get 


dtipl(t;uj,N) = 


r*+oo 


dtip 2 (t-,u,N) 


r»+oo 


yphjj 


d(C N {Ln (C)) 2 exp(-C). 


But 


F?(z) = - 


\/ 2 cj , 


where 


N+2n 


*) = - E ^Ww)') 


r+oc -I 

d(C N {Ln (0) exp(-C) = -^=exp {-z)Q*(z), 


d? 


3=0 


(5.2) 


(5.3) 


is a polynomial of degree N + 2n. 

In contrast with what happened in the TDPT case, to obtain a rational extension we need 
now to eliminate the exponential factor, which necessitates fixing Ai = — ||-0 n || 2 . Equation (3.6) 
then becomes 

(x; w, N, - UJ 2 ) = V(x-u,N)-2[ log (Ef (*))]" 

i / 2 d [z N+1 / 2 (L 1 //(z)) 2 \ 

= V{x-,u,N)-4uz —y J, 

where V^ n 2 \x; uj, N, — ||V ; n|[ 2 ) constitutes a rational extension of V(x\uj,N). Since (z) is a 
negative, strictly increasing function of z, such that 


lim Ff (z) = 0, 

x-H-oo 
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Qn (z) keeps a constant strictly negative sign, i.e., Q% ( 2 ) has no zero on the positive half-line. 
This shows that V^ n \x;uj,N, — \\il > n || 2 ) is also regular. 

Let us note too that [18, 50] 


or 


?N 


Qn( 0 ) = - 


1 f + °° o 

( 0 ) = —7= / d(C N (L%( 0 ) exp(-C) = - 

\/ 2 w jo 

(n + IV)! 


1 (n + lV)! 

a/ 2 cJ n\ 


< 0 , 


n\ 


The eigenfunctions of V^ n2 \x;uj,N, — ||^ n || 2 ) are given by (see equation (3.7)) 

\x-,u,N, -HV’nll 2 ) = [E n (u) - E'fc(w)]V , fc(x; to, N) - W(ip n ,4k \ x ) — 

Ti (x; cj, N, — ||V’n|| 2 ) 


where 


^(x;u,N,-\ 


Jo ^n(*i jV ) ~ 

ip n (x;u,N) 


F = _J_ -(JV+l/2)/2 /_ 

Mw,N) ^ eXPl * A ] L*{z) 


and 


W(4> n ,4 k | x) = v/2w^ JV+1 exp(-2;)L^ A ,(2), 


with 


< fc (z) = L™(z)L%(z) - L%(z)L%+\z). 
Consequently 

V’jJ \x;u,N,-\\4> n \\ 2 ) = [E n {uj) — Ek(u])\4>k(x]U), N) 


- 2 ^.( 3 *+W exp(—z/ 2 ) w , 


that is, 


_7X** a ) /_Ar 11 .,. n2\ _ 2wz( Ar+1/2)/2 exp(-2;/2) 7 ( ri 2)^^ 

(^)W, IV, HV'nll ) Q N (z) F k,N (^0’ 


where 


4jv(~) = ( n - fe ) L fc (*)Qn (*) - zJV+li n,t( z ) i n (*)■ 

V>), n (a:; w, IV, — HV’nll 2 ) decreases exponentially to zero at infinity and also tends to zero at 
the origin (see equation (5.5)). It is therefore an admissible eigenstate of V^ n \x\oj, N, — ||V , n|| 2 ) 
for every k 7 ^ n. 

Note that in this case 

4f\x-u,N,-W4nW 2 ) = lM n )(s; W ,JV,-||^„|| 2 ) OC z (^v+i/ 2)/2 eM z/ 2 )^l (5.7) 

Vn \ z ) 


(5.4) 


(5.5) 


k 7 ^ n, 


(5.6) 
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is not an eigenstate since it is not normalizable. This implies that V^ n2 \x] u>, N, — ||^> n || 2 ) and 
V(x;u,N) are only quasi-isospectral, the two-step confluent DBT A(Wv) o A(w v ) being now 
state-deleting. 

The orthogonality conditions between eigenstates imply that the set of (z) with k / n, 
constitute a family of orthogonal polynomials, indexed by k E N, on the positive half-line with 
respect to the mesure 



z N exp(— z) 

W (*)) 2 


5.2 Examples 
5.2.1 The n = 0 case 

The n = 0 case gives no new extended potentials with respect to previously known results. 
Indeed, the corresponding confluent two-step DBT is obtained by applying successively the one- 
step DBT y4(tco) and A(ITg 0 ^). The first one is the usual SUSY partnership and, due to the 
shape invariance property of the isotonic potential, we have 

U (0) (x; u, N ) = V(x; u, N + 1) + 2w, 

where admits the following eigenstates 

uj, N) ~ ipk-i(x; w,N + 1), 

for the respective energies Ek( w), k > 1. 

The second DBT A(W^ 0) ) is associated to the seed function T ( 0 0) (x;u;,N,-||^o|| 2 ), which 
is a formal eigenfunction of (x; &, N) for the eigenvalue Eo = 0 and which is then in the 
disconjugacy sector of V^(x; u, N). But 


t[, 0 ) (x;u;, JV, Ai) = ip§\x;w,N) ( Ai + [ dt —~— --^ ) , 

V Jx ° {4\t^N)) 2 j 

where 

ipQ 0 \x; U), N) ~ l/ipo(x\u,N) = z -^+i/ 2 )/ 2 e +z /2 _ 

By taking xo = 0 and Ai = — 11 -0011 3 , we obtain 

Ti 0) (x;u;,iV,-||^o|| 2 ) = exp(-z/2)Q?(z), (5.8) 

V zcj 

where 

N N l 

3=0 1=0 

which gives 

<P< 0 ) (x; W ,JV,-||V.„|| 2 ) = —Y ex p(- 2 /2) z -(»+>/2)/2^f!, (5.9) 

l=Q l. 
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Equation (5.9) shows that Tq 0 ^(x; uj, N, — HV’oll 2 ) is a formal eigenfunction of V^(x;u,N) = 
V (x; uj, N+1)+2uj that tends to zero at infinity and diverges at the origin. Such an eigenfunction 
is unique (up to a constant multiplicative factor) at a given eigenvalue (here Eq = 0) [ 6 , 30]. 

Consequently ^^(x; uj, N, — ||V’o|| 2 ) coincides necessarily (up to a constant multiplicative 
factor) with a type II seed function [22] of V (x; uj, N + 1) + 2cj, which is a formal eigenfunction 
for the eigenvalue E— \ + 2uj, namely 

^|, 0 ) (x;a;,iV,-||^o|| 2 ) ~ <t>N ,-(x; uj, N + 1 ) = z~ l ' N+1/2)/2 e~ z/2 LJ f N ~ l (z), 

where (j} n - is defined by 

<j> n -(x\uj,a ) = i/j n (x;uj, -a). 


Note that this result can be recovered directly by using the identity (see equation (5.1)) 

N i 

D A 'Ejr. 

1=0 

in equation (5.9). 

Applying then the DBT A(Wq°^) = A(vn,-(x;uj,N + 1)) (vjy,- = —(j) , N _/(j)N,-) to V^\x\ 
uj,N) = V(x',uj, N + 1) + 2uj, we arrive at a two-step confluent extended potential 

V ( ° 2) (x; u, N, -HV’oll 2 ) = V^ N ’~\x-uj, IV + 1) + 2ca, 

which coincides (up to a constant shift) with the usual extensions obtained with type II seed 
functions [ 22 ], 

Explicitly, we get 


p(° 2 )(x;u;,IV,-||V’o|| 2 ) = V(x;uj,N) + ^^ 1/2 |- 


N 


N+l/2 


E Ei'- 


1=0 


and 


^° 2) (x; 

uj,N, 

-liv’d! 2 ) 


z (iV+i/2)/2 exp( 

; AT 






£ d//! 






1=0 

In the IV = 

= 1 case, this gives 



P(° 2 )(x: 

;w,l,- 

-IIV’oll 2 ) 

a; 2 2 

= T x + 

3 

4x 2 

4uj 

cux 2 + 2 


> 1. 


16cu 


= V^ 1, \x;u, 2) + 2 uj 


and 


1° } (x; uj, 1, — 


2 \ 2wz 3 / 4 exp(—z/2) ~(o 2 ) ( ^ 

) ~ _ , i L k. 1 (Ah 


z + 1 


k > 1, 


“ kLl(z)(z + 1) + z 2 Ll_ 1 (z). 


with 
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5.2.2 The n = 1 case 

We have in this case 


N 


Qi{z) = - z N+ 2 + Nz n+1 — (N + 1 )! ^2 


j =o 


and 


AT,—HVhH 2 ) = V(x;v,N)~ 4o T =) [(N + l)(N + 1 ) - (N + \)z] 

vi \ z ) 


— 4a; 


+ 1 - z) s 


N-l 


[Qn4 2 

The eigenstates of are (k / 1) 


(N + 2 )z N+1 - N(N + 1 )z N + {N + 1)! ^ — | . 

i=o 


riV / \r./v/ \ 

4 } (aj; Cc;, iV, — ||V >1 II 2 ) = (1 - k)^k(x; u, N ) - Z ( 3N + 5 / 2 )/ 2 exp(-z/2) 1,A nAf 1 - , 

Vl GJ 


that is, 


4 ) (x;w,iV,-||V’i|| 2 ) = 


2 w 2:( Ar+1//2 )/ 2 exp(— z/2) ~(i 2 ) 


Qf(*) 




where 


4$(*) = ^ +1 (JV + 1 - z) a i£Y(*) 


2 riV+1/ 


AT 


+ ( A :^ 2 - (IVfc + l )^ 1 + (fc - l)(iV + 1)! ^ — ) L%(z). 

3=0 J ' 


In particular, for N = 1 

Qj(z) = — z 3 + z 2 — 2z — 2 


and 


V {l2) (x-u, l,-||V’i|| 2 ) = V(x;u, 1) + 6 w ^3 


z 2 -4 
z 2 + 2z + 2 


— 20 a; 


5z 2 - 4z - 2 

(z 3 - z 2 + 2z + 2) 2 


+ 2 cj, 


or 


y (l 2 ) (x;w,l,-||^i|| 2 ) 


a; 2 o 3 a; 2 x 4 — 16 

-'T -I- - -I- 1 - 

4 4x 2 *" w 3 x 6 — 2 cj 2 x 4 + 8 ujx 2 + 16 

5w 2 x 4 — 8cjx 2 — 8 

— 320a;-2- 

(a; 3 x 6 — 2 u; 2 £ 4 + 8o;.t 2 + 16)" 


Moreover (k ^ 1) 



2ujz 3 ^ exp(—z/ 2 ) ~(i 2 ), . 
z 3 ~z 2 + 2z + 2 fc ’ 4 [Zh 


where 


L^i(z) = [kz 3 - (fc + l)z 2 + 2(k - 1 )z + 2{k - 1 )\L\(z) + z 2 (2 - z) 2 L^_ 1 (z). 



Confluent Chains of DBT: Enlarged Shape Invariance and New Orthogonal Polynomials 21 


The polynomial 

L^\z) = -(z 2 + 2z + 2) 

is associated to the ground state of 0 , 1 , — ||?/h|| 2 ) and 

L 2 I ) (z) = - \ (- 4 + 4^ 2 - 12 ), 

which has a single zero at \/2 on ] 0 , +oo[, corresponds to the first excited eigenstate. 


5.3 Shape invariance of the two-step confluent rational extensions 
of the isotonic potential 

V(x',uj,N) is a translationally shape invariant potential with a SUSY partner 
Vsvsy(x- u, N ) = V<® (x; u,N) = U(x; u, N + 1) + E^u) 


and (the coefficient is readily established by a direct calculation) 

^n\x]ui,N) = = -\ / 2wV'n-i(x;w,iV + 1). 

V’o (x;w,IVj 


(5.10) 


Considering the SUSY partner of (x; u, N, — HV’nll 2 ), let us deal separately with the 
n > 1 and n = 0 cases. For n > 1, this partner is given by 


^SUSyO^ HIlM 2 ) = U^ n2) (x; w, N , —HV’nll 2 ) - 2( log 4^ ] (x; w, N, - ||^ n || 2 ))" (5.11) 




= V(x;u,N) - 2[log fa n (x]u,N)^\x;u,N,-\\^ n \\ 2 )^ n \x;u,N, — 11^11 2 ))]", 


where 


V>o n ) (x;u;,!V,-||V’n|| 2 ) = E n (u)ip 0 {.x;u;, N) - 


W (V’n, V’O | z) 


(x;U,N, -HV'nll 2 ) 


(5.12) 


Inserting (5.12) into (5.11), we get 

^susY^^’^-IIV’nll 2 ) = V (x; u, N) - 2(log'ip 0 (x-uj,N)) r 
- 2 | log 
= V(x; u, N + 1) + Ei (ui) 

{log (2uuj^^(x-,uj, N, -HV’nll 2 ) - V2ui/j n -i{x',u, N + 1)^ ip n (x; u , N) | . 


E n (o;)^)(x;u;,IV,-||^n|| 2 ) + ) ij) n (x\ 10 , N) 

ipo{xiu,N) 


-2 


More precisely, on using (5.6), 

^susy(x; w , 2V, — 11 -0^ 11 2 ) = V(x;uj,N + 1) — 2(log'0n-i(x; w, N + 1))" + E 1 (u) 


~ 2 < log 


- 


= V {n - 1] (x; w, N + 1) + E 1 (w) 


, N V’n(x;w,IV)V’n-i(x;a;,IV +1) 




l on 


- 2 ^ log 


jv Vn(x;o;, N)il> n _i(x-,u, N + 1) 


- 


y/2uj’ 


um 


/lp n -l(x-,UJ,N + 1) 


/lp n -l(x-,UJ,N + 1) 
















22 


Y. Grandati and C. Quesne 


If, for an appropriate constant C, the following condition 


f n^ _ i/’n(x;u},N)i/j n - 1 (x-,u,N + l) = cf n+i^ 




um 


is satisfied, then 


vN, n ip n (x;uj,N)ip n _ i(x;w,IV + l)\ / 

F n (z) - 1 = - /'lp n - 1 (x-,UJ,N + 1) 

\/2wn J ' 

= C*%r i 1) (x-,w,N + l,-||^„_i|| 2 ) 
and we obtain an enlarged shape invariance property 

V^ Y {x;u;,N,-U n \\ 2 ) = (x; a;, iV + 1 ,-||^n-if) + 

The preceding condition (5.13) can be rewritten as 
A(z) — B(z) = 0, VzG]—1,1 [, 

with 

-y r+°° 


(5.13) 


(5.14) 

(5.15) 


A(z) = F?(z) - CF£_\\z) = ^=l dc(cC N+1 {L»+\ C )) 2 - C N K( C)) 2 ) exp(-C) 


and 


£(*) = 


tpnix] w, N)ip n -i(x] cu,N + 1) z N+l e- z L%{z)L%+l{z) 


\/2 ujn y/2um 

Equation (5.15) is equivalent to the set of conditions 

%-A (z) = | -B(z ) 


(5.16) 


and 


A(0) = 5(0). (5.17) 

Here F A(z) is given by 

^^Z A (z) = (( L n ( z )f “ Cz { L n-l(z)) 2 )z N exp (- Z ). 

As for f B(z), using the derivation formula [18, 50] 

4i«( z ) = -iy«( z ), 

it can be expressed as 

V^un-^Biz) = -z N exp (~z)L%(z) (z^L%(z) + (N + 1 — z)~L^(z)J 

-z N+1 eM-z){L^l(z)) 2 . (5.18) 

The differential equation satisfied by the Laguerre polynomials is [18, 50] 
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which, inserted in equation (5.18), yields 

V / 2 u-^B(z) = {(L»(z)) 2 - ^z{L»+\z)) 2 )z N eM-z). 
Hence, to satisfy the first condition (5.16), we must choose 



With such a choice for C, the second condition (5.17) is then also fulfilled because it is 
obvious that B{ 0) = 0 and, from equation (5.4), it results that 

H(0) = J F n Ar (0)-^_+ 1 (0) = 0. 

We conclude that for any n > 1, the two-step confluent extension VguSY^iw,IV, ~ HV’nll 2 ) 
satisfies the enlarged shape invariant property (5.14). 

Such a property is in particular satisfied if n = N = 1, as it can be directly verified. In this 
case, we indeed have 

L { o\z) = 1, L^\z) = 2-z, 


and C = 1, so that the identity (see equations (5.2) and (5.3)) 


A(z) = Fl(z)-FZ(z) = 


y/2u 


ex' 


v{~z){Q\{z) - Ql{z)) 


-4= exp(^) | v ig ( 3 2 ) -T.F (^<2 - zf) | = «,-* = B(z 

v 2j l 4?, dz> ' dzl > I 


F( z *)-<r d ’ 

^ h 

is fulfilled, which implies that 

^susy(‘ t ; w ’ “IIV'ill 2 ) = W ((0)2) (a^;c^,2, — IIV^oll 2 ) +^i(w). 

Let us now turn ourselves to the n = 0 case. As noted above (see equation (5.7)), the formal 
eigenfunction 


tf" ) (-'CW,A^-||V’o|| 2 ) = 


1 


^ 0) (x;a;,lV,-||^o|| 2 ) 


associated to the eigenvalue Eq = 0, is not an eigenstate of (x; ui, N, — ||V’o|| 2 ) and the 
ground state of this potential is actually given by (x; u, N, -||^o|| 2 ) • The SUSY partner 


therefore reads 
do 2 ) 


V s\jsy{ x 'i u ’ N ’- IIV’oll 2 ) = U ( ° 2) (x;a;, JV, — ||^o|| 2 ) -2(log^° } (x;w, N, — H^o II 2 )) ^ 


= V(x;u,N) -2 log (x, N)^ ] (x- ui, N, - } (x; u, N, -||V’o|| 2 )) 


2U7(o 2 b 


4° ) (-'c;w, A r ,-||V’o|| 2 ) = -Ei(w)V’i(x;o;,lV) - 


W(ip 0 ,^ 1 1 x) 


Tq 0) [x] uj, n, — liv’d! 2 ) 


where 
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Using the shape invariance property of V, we can then write 

®y(*;w»^»-H^)II 2 ) = V{x ] u,N + l) + E 1 {u) 

log (—20/01 (x; u, A r )4'g°' ) (x; u, N, -||0o|| 2 ) ~W(^ ’0,01 [x)) 


- 2 


that is, with equation (5.8) and (5.10), 

^s ( usy(*; w »^-H^|| 2 ) = V{x-,U,N+1 ) + E 1 (uj) - 2[log0 o (^,iV)]" 
log {L?(z)z~( n+1 W 2 exp(—z/2)Q^ r (z) - 0 o (x; u, N + 1)) 


- 2 


= V {0 \x]uj,N + 1) + E 1 (uj ) 

- 21 log [.-^+ 3 / 2 )/ 2 eX p(- z /2)(Lf (z)Q$(z) 


— z 


N +1 


If there existed some constant C such that 

L?(z)Q$(z) - z N+1 = CQq +1 (z), (5.19) 

then we would get 

^SUSyM, = V^ 2 \x-u,N + 1) + E 1 (u). 

However, it can be readily seen that equation (5.19) cannot be satisfied by any C, so that 
we conclude that we do not have any strict nor enlarged shape invariance for the confluent 
extension TA° \ 


6 Final comments 

Using two-step confluent chains of DBT, we have generated new families of orthogonal poly¬ 
nomials, associated to novel regular rational extensions of the isotonic and TDPT potentials, 
exhibiting an enlarged shape invariance property. Interestingly, in the second case, the ortho¬ 
gonal polynomials depend on a free parameter that can be modulated continuously, a feature 
already encountered for other extensions based on para-Jacobi polynomials. 

Considering chains of arbitrary order and the possibility of obtaining more general families 
of orthogonal polynomials subjected to multi-parameter dependence would be a very interesting 
topic for future investigation. 

As a final point, it is worth observing that apart from the Darboux (or differential) approach 
to the construction of (quasi-)isospectral families of potentials, considered in the present paper, 
there also exist the Abraham-Moses [1] and Pursey [40] (or integral) approaches to the same, 
coming from the inverse scattering Gel’fand-Levitan or Marchenko technique. Relations between 
both types of methods have been extensively studied by several authors (see, e.g., [5, 33, 38, 46]), 
showing that they sometimes lead to the same results, but are in general inequivalent. The 
differential procedure being more convenient and easier to use in quantum mechanics has been 
preferred here. Studying the relation between our results and those of the Abraham-Moses 
method would, however, be an interesting open question for future work. 


References 

[1] Abraham P.B., Moses H.E., Changes in potentials due to changes in the point spectrum: anharmonic 
oscillators with exact solutions, Phys. Rev. A 22 (1980), 1333-1340. 







Confluent Chains of DBT: Enlarged Shape Invariance and New Orthogonal Polynomials 25 


[2] Adler M., Moser J., On a class of polynomials connected with the Korteweg-de Vries equation, Comm. 
Math. Phys. 61 (1978), 1-30. 

[3] Adler V.E., On a modification of Crum’s method, Theoret. and Math. Phys. 101 (1994), 1381 1386. 

[4] Bagchi B., Grandati Y., Quesne C., Rational extensions of the trigonometric Darboux-Poschl Teller poten¬ 
tial based on para-Jacobi polynomials, J. Math. Phys. 56 (2015), 062103, 11 pages, arXiv: 1411.7857. 

[5] Baye D., Phase-equivalent potentials for arbitrary modifications of the bound spectrum, Phys. Rev. A 48 
(1993), 2040-2047. 

[6] Berezin F.A., Shubin M.A., The Schrodinger equation, Kluwer, Dordrecht, 1991. 

[7] Bermudez D., Fernandez C. D.J., Fernandez-Garcia N., Wronskian differential formula for confluent super- 
symmetric quantum mechanics, Phys. Lett. A 376 (2012), 692-696, arXiv:1109.0079. 

[8] Burchnall J.L., Chaundy T.W., A set of differential equations which can be solved by polynomials, Proc. 
London Math. Soc. S2-30 (1930), 401-414. 

[9] Calogero F., Yi G., Can the general solution of the second-order ODE characterizing Jacobi polynomials be 
polynomial ?, J. Phys. A: Math. Theor. 45 (2012), 095206, 4 pages. 

[10] Carinena J.F., Ramos A., Integrability of the Riccati equation from a group-theoretical viewpoint, Inter- 
nat. J. Modern Phys. A 14 (1999), 1935-1951, math-ph/9810005. 

[11] Carinena J.F., Ramos A., Fernandez C. D.J., Group theoretical approach to the intertwined Hamiltonians, 
Ann. Physics 292 (2001), 42-66, math-ph/0311029. 

[12] Contreras-Astorga A., Fernandez C. D.J., Supersymmetric partners of the trigonometric Poschl-Teller po¬ 
tentials, J. Phys. A: Math. Theor. 41 (2008), 475303, 18 pages, arXiv:0809.2760. 

[13] Cooper F., Khare A., Sukhatme U., Supersymmetry in quantum mechanics, World Scientific Publishing 
Co., Inc., River Edge, NJ, 2001. 

[14] Crum M.M., Associated Sturm-Liouville systems, Quart. J. Math. Oxford Ser. (2) 6 (1955), 121-127, 
physics/9908019. 

[15] Darboux G., Sur une proposition relative aux equations lineaires, C. R. Acad. Sci. Paris 94 (1882), 1456- 
1459. 

[16] Darboux G., Legons sur la theorie generale des surfaces et les applications geometriques du calcul 
infinitesimal. II, 2nd ed., Gauthier-Villars, Paris, 1915. 

[17] Duran A.J., Perez M., Admissibility condition for exceptional Laguerre polynomials, J. Math. Anal. Appl. 
424 (2015), 1042-1053, arXiv:1409.4901. 

[18] Erdelyi A., Magnus W., Oberhettinger F., Tricomi F.G., Higher transcendental functions, McGraw-Hill, 
New York, 1953. 

[19] Fernandez C. D.J., Salinas-Hernandez E., The confluent algorithm in second-order supersymmetric quantum 
mechanics, J. Phys. A: Math. Gen. 36 (2003), 2537-2543, quant-ph/0303123. 

[20] Fernandez C. D.J., Salinas-Hernandez E., Wronskian formula for confluent second-order supersymmetric 
quantum mechanics, Phys. Lett. A 338 (2005), 13-18, quant-ph/0502147. 

[21] Fernandez C. D.J., Salinas-Hernandez E., Hyperconfluent third-order supersymmetric quantum mechanics, 
J. Phys. A: Math. Theor. 44 (2011), 365302, 11 pages, arXiv:1105.2333. 

[22] Gomez-Ullate D., Grandati Y., Milson R., Extended Krein-Adler theorem for the translationally shape 
invariant potentials, J. Math. Phys. 55 (2014), 043510, 30 pages, arXiv: 1309.3756. 

[23] Gomez-Ullate D., Kamran N., Milson R., An extended class of orthogonal polynomials defined by a Sturm- 
Liouville problem, J. Math. Anal. Appl. 359 (2009), 352-367, arXiv:0807.3939. 

[241 Grandati Y., Solvable rational extensions of the isotonic oscillator, Ann. Physics 326 (2011), 2074-2090, 
arXiv:1101.0055. 

[25] Grandati Y., Multistep DBT and regular rational extensions of the isotonic oscillator, Ann. Physics 327 
(2012), 2411-2431, arXiv: 1108.4503. 

[26] Grandati Y., New rational extensions of solvable potentials with finite bound state spectrum, Phys. Lett. A 
376 (2012), 2866-2872, arXiv: 1203.4149. 

[27] Grandati Y., A short proof of the Gaillard-Matveev theorem based on shape invariance arguments, Phys. 
Lett. A 378 (2014), 1755-1759, arXiv:1211.2392. 

[28] Grandati Y., Berard A., Rational solutions for the Riccati-Schrodinger equations associated to translation- 
ally shape invariant potentials, Ann. Physics 325 (2010), 1235-1259, arXiv:0910.4810. 



26 


Y. Grandati and C. Quesne 


[29] Grandati Y., Berard A., Comments on the generalized SUSY QM partnership for Darboux-Poschl-Teller 
potential and exceptional Jacobi polynomials, J. Engrg. Math. 82 (2013), 161-171. 

[30] Hartman P., Ordinary differential equations, John Wiley & Sons, Inc., New York London - Sydney, 1964. 

[31] Keung W.-Y., Sukhatme U.P., Wang Q.M., Imbo T.D., Families of strictly isospectral potentials, J. Phys. A: 
Math. Gen. 22 (1989), L987-L992. 

[32] Krein M.G., On a continual analogue of a Christoffel formula from the theory of orthogonal polynomials, 
Dokl. Akad. Nauk SSSR 113 (1957), 970-973. 

[33] Luban M., Pursey D.L., New Schrodinger equations for old: inequivalence of the Darboux and Abraham- 
Moses constructions, Phys. Rev. D 33 (1986), 431-436. 

[34] Matveev V.B., Generalized Wronskian formula for solutions of the KdV equations: first applications, Phys. 
Lett. A 166 (1992), 205-208. 

[35] Matveev V.B., Positons: slowly decreasing analogues of solitons, Theoret. and Math. Phys. 131 (2002), 
483-497. 

[36] Messiah A., Mecanique quantique, Vol. 1, Dunod, Paris, 1959. 

[37] Mielnik B., Nieto L.M., Rosas-Ortiz O., The finite difference algorithm for higher order supersymmetry, 
Phys. Lett. A 269 (2000), 70-78, quant-ph/0004024. 

[38] Nieto M.M., Relationship between supersymmetry and the inverse method in quantum mechanics, Phys. 
Lett. B 145 (1984), 208-210. 

[39] Odake S., Sasaki R., Infinitely many shape invariant potentials and new orthogonal polynomials, Phys. 
Lett. B 679 (2009), 414-417, arXiv:0906.0142. 

[40] Pursey D.L., New families of isospectral Hamiltonians, Phys. Rev. D 33 (1986), 1048-1055. 

[41] Quesne C., Exceptional orthogonal polynomials, exactly solvable potentials and supersymmetry, J. Phys. A: 
Math. Theor. 41 (2008), 392001, 6 pages, arXiv:0807.4087. 

[42] Quesne C., Solvable rational potentials and exceptional orthogonal polynomials in supersymmetric quantum 
mechanics, SIGMA 5 (2009), 084, 24 pages, arXiv:0906.2331. 

[43] Quesne C., Higher-order SUSY, exactly solvable potentials, and exceptional orthogonal polynomials, Modem 
Phys. Lett. A 26 (2011), 1843-1852, arXiv: 1106.1990. 

[44] Quesne C., Novel enlarged shape invariance property and exactly solvable rational extensions of the Rosen- 
Morse II and Eckart potentials, SIGMA 8 (2012), 080, 19 pages, arXiv:1208.6165. 

[45] Quesne C., Revisiting (quasi-)exactly solvable rational extensions of the Morse potential, Internat. J. Modem 
Phys. A 27 (2012), 1250073, 18 pages, arXiv:1203.1811. 

[46] Samsonov B.F., On the equivalence of the integral and the differential exact solution generation methods 
for the one-dimensional Schrodinger equation, J. Phys. A: Math. Gen. 28 (1995), 6989-6998. 

[47] Samsonov B.F., New possibilities for supersymmetry breakdown in quantum mechanics and second-order 
irreducible Darboux transformations, Phys. Lett. A 263 (1999), 274-280, quant-ph/9904009. 

[48] Schulze-Halberg A., Wronskian representation for confluent supersymmetric transformation chains of arbi¬ 
trary order, Eur. Phys. J. Plus 128 (2013), 69, 17 pages. 

[491 Sparenberg J.-M., Baye D., Supersymmetric transformations of real potentials on the line, J. Phys. A: Math. 
Gen. 28 (1995), 5079-5095. 

[50] Szego G., Orthogonal polynomials, American Mathematical Society Colloquium Publications, Vol. 23, 4th ed., 
Amer. Math. Soc., Providence, R.I., 1975. 

[51] Vein R., Dale P., Determinants and their applications in mathematical physics, Applied Mathematical Scien¬ 
ces, Vol. 134, Springer-Verlag, New York, 1999. 



